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ABSTRACT 

It  is  well  known  that  the  growth  rate  of  instability  waves  of  a  two-dimensional  free 
shear  layer  is  reduced  greatly  at  supersonic  convective  Mach  numbers.  In  previous  works,  it 
has  been  shown  that  new  wave  modes  exist  when  the  shear  layers  are  bounded  by  a  channel 
due  to  the  coupling  effect  between  the  acoustic  wave  modes  and  the  motion  of  the  mixing 
layer.  The  present  work  studies  the  simultaneous  propagation  of  multiple  stability  waves 
using  numerical  simulation.  It  is  shown  here  that  the  co-existence  of  two  wave  modes  in 
the  flow  field  can  lead  to  an  oscillatory  growth  of  disturbance  energy  with  each  individual 
wave  mode  propagating  linearly.  This  is  particularly  important  when  the  growth  rates  of 
the  unstable  waves  are  small.  It  is  also  shown  here  that  the  propagation  of  two  neutrally 
stable  wave  modes  can  lead  to  a  stationary  periodic  structure  of  r.m.s  fluctuations.  In 
the  numerical  simulations  presented  here  the  forced  wave  modes  are  propagating  at  same 
frequency  but  with  different  phase  velocities.  In  order  to  track  the  growth  of  each  wave 
mode  as  it  propagates  downstream,  a  numerical  method  which  can  effectively  detect  and 
separate  the  contribution  of  individual  wave  is  given.  It  is  demonstrated  that  by  a  least 
square  fitting  of  the  disturbance  field  with  eigenfunctions  the  amplitude  of  each  wave  mode 
can  be  found.  Satisfactory  results  as  compared  to  linear  theory  are  obtained. 
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under  NASA  Contract  NASl-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer 
Applications  in  Science  and  Engineering,  NASA  Langley  Research  Center.  Hampton,  VA  23665, 
USA. 
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I.  INTRODUCTION 


In  this  paper  we  present  results  of  numerical  simulations  of  a  confined  supersonic 
mixing  layer  (Figure  1).  It  is  well  known  that  at  high  supersonic  convective  Mach  numbers 
(Mach  number  in  the  reference  frame  of  the  large  scale  structure)  the  familiar  Kelvin- 
Helmholtz  instability  wave  of  a  free  mixing  layer  becomes  stabilized  due  to  the  increased 
compressibility  of  the  flow.  It  was  found  in  recent  experiments^’^’*  as  well  as  in  numerical 
simulations'*’®’®  that  the  spreading  rate  of  a  mixing  layer  at  a  supersonic  convective  Mach 
number  is  a  factor  of  four  or  five  smaller  than  at  a  low  subsonic  convective  Mach  number. 
As  a  result,  large  scale  roll-ups  of  the  vorticity  layer,  a  standard  feature  of  low  speed  free 
shear  layers,  are  absent  and  the  mixing  of  the  two  streams  is  inefficient.  In  efforts  to  explain 
these  observations,  a  number  of  recent  studies  on  the  stability  of  the  free  shear  layer  using 
inviscid  quasi-parallel  linear  theory  have  been  carried  out^’®.  These  studies  showed  that 
the  intrinsic  instability  wave  of  the  shear  layers  has  a  much  rediiced  growth  rate  when  the 
convective  Mach  number  becomes  supersonic.  It  was  suggested  that  this  reduction  of  the 
unstable  wave  growth  is  directly  responsible  for  the  small  spreading  rate  of  the  supersonic 
mixing  layers  observed  in  the  experiments  and  direct  numerical  simulations. 

At  the  same  time,  it  was  also  found  that  a  bounded  free  shear  layer  behaves  differently 
from  an  unbounded  one  at  high  supersonic  convective  Mach  numbers^’’®’".  These  studies 
revealed  that  the  same  shear  layer,  when  bounded  by  a  channel,  consisting  of  ffat  walls  at 
the  top  and  bottom,  supports  new  instabilities  due  to  the  coupling  between  the  acoustic 
wave  modes  of  the  channel  and  the  motion  of  the  shear  layer.  A  thorough  treatment 
of  normal  modes  associated  with  a  bounded  shear  layer  was  given  by  Tam  and  Hu®. 
Systematic  calculations  of  normal  mode  solutions  showed  that  four  families  of  waves  exist. 
In  their  naming  convention,  Class  A  and  Class  B  modes  are  the  unstable  waves,  Class  A 
waves  are  related  to  the  wave  reflections  from  the  lower  wall  and  Class  B  waves  are  related 
to  the  reflections  from  the  upper  wall.  In  addition,  there  are  two  families  of  neutrally 
stable  waves.  They  were  called  Class  C  and  Class  D  acoustic  waves  since  they  are  related 
to  the  acoustic  reflection  off  both  walls.  Typical  dispersion  relations  of  the  four  waves  and 
the  growth  rates  of  A  and  B  waves  are  given  in  Figure  2  where  the  calculations  have  been 
made  for  two-dimensional  waves  at  Mj  (Mach  number  of  the  upper  stream)=  3.5.  M2 
(Mach  number  of  the  lower  stream)=  1.2  and  the  sound  speeds  ratio  of  the  two  streams 
01/02  =  1.2.  In  Figure  2  a  subscript  has  been  used  to  indicate  the  corre.sponding  mode 
number.  It  is  to  be  noted  that  although  the  D  waves  have  negative  phase  velocities,  all 
waves  have  positive  group  velocities  and  hence  are  downstream  propagating  waves. 

’When  disturbances  are  initiated  upstream,  it  is  expected  that  the  wave  with  the 
largest  growth  rate  dominates  downstream.  Observation  of  the  spatial  growth  rates  of 
the  A  and  B  waves  indicate  that  at  certain  frequencies  the  A  waves  are  donunant  while 
at  other  frequencies  the  B  waves  are  the  dominant  waves.  Yet  the  A  and  B  waves  have 
comparable  growth  rates.  In  numerical  simulations  unless  the  upstream  perturbati'^n  is  an 
eigenfunction  of  a  single  wave  mode  it  is  iikeiy  that  both  the  A  and  B  unstable  waves  as 
well  as  the  neutrally  stable  waves  are  excited  and  propagate  downstream.  Consequently, 
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simulations  of  a  bounded  mixing  layer  could  be  significantly  different  from  the  unbounded 
counterpart  in  two  respects.  First,  since  the  A  and  B  waves  have  comparable  growth  rate, 
both  waves  are  important  in  downstream  propagation.  Second,  becavise  the  growth  rates 
of  the  \mstable  A  and  B  waves  are  generally  small  at  high  convective  Macli  numbers,  no 
large  scale  vorticity  roll-ups  are  present  and,  as  a  result,  the  excited  neutrally  stalde  waves 
may  also  constitute  an  important  part  in  the  downstream  flow  field. 

The  present  work  is  concerned  with  the  cases  when  multiple  wav'es  are  forced  at  the 
same  frequency.  One  objective  of  this  paper  is  to  show  that  the  co-existence  of  multiple 
propagating  waves  can  lead  to  an  appearance  of  oscillatory  growth  of  the  disturbance 
flow  field  even  though  the  waves  are  growing  linearly.  A  second  objective  of  this  paper  is 
to  construct  a  numerical  method  that  can  detect  the  growth  of  individual  wave  a.s  they 
propagate. 

In  many  experimental  and  computational  studies  it  is  essential  to  know  the  growtli  or 
decay  rate  of  the  propagating  waves.  For  example,  in  the  studies  of  turbulent  transitions  of 
plane  mixing  layers’^’’^  ,  the  spatial  evolutions  of  the  fundamental  and  subharmonic  waves 
are  monitored  as  the  waves  propagate  downstream.  The  interplay  of  the  fundamental  and 
subharmonic  waves  is  studied  by  their  growth  and  saturation  as  they  propagate.  Usually 
the  waves  are  identified  by  their  frequencies  using  time  series  spectrum  analysis  of  the 
disturbances  at  chosen  downstream  locations.  The  assumption  is  that,  for  plane  mixing 
layers,  the  only  dominant  growing  (unstable)  wave  mode  is  the  Kelvin-Helmholtz  wave. 
Thus  the  frequency  spectrum  naturally  indicates  the  strength  of  the  K-H  wave  at  any  given 
frequency.  For  this  reason,  the  local  wavenumber  of  the  disturbances  can  also  be  estimated 
from  the  phase  difference  of  simultaneous  measurements  of  the  time  series  at  two  nearby 
points  in  space' In  this  way  the  dispersion  relationship  can  be  established  experimentally. 
However,  when  multiple  wave  modes  of  comparable  strength  are  propagating  at  the  same 
frequency  but  with  different  phase  velocities,  the  growth  or  decay  of  each  individual  wave 
mode  can  not  be  distinguished  from  an  examination  of  the  frequency  spectrum  alone.  One 
method  of  tracking  each  wave  is  to  use  a  cross-bispectrum  analysis'^.  In  this  paper  a 
method  based  on  the  eigenfunction  expansion  is  proposed.  It  is  shown  numerically  that  by 
an  expansion  of  the  disturbances  in  the  form  of  eigenfunctions  it  is  possible  to  track  and 
separate  the  growth  of  each  individual  mode. 

In  section  2  we  give  the  formulation  of  the  problem  and  a  description  of  the  numerical 
methods.  The  method  of  eigenfunction  expansion  for  detecting  wave  ami)litu<les  is  also 
described.  In  section  3  results  of  numerical  simulations  and  eigenfunction  expansions  an’ 
presented.  Two  ca.ses  are  considered  and  discu.ssed.  Section  4  contains  the  conclusions. 


II.  FORMULATION 

A.  Governing  equations  and  numerical  .simulation 

We  consider  a  two-dimensional  mixing  layer  formed  between  two  streams  of  uniform 
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parallel  flow  confined  inside  a  channel  (Figure  1).  The  coordinate  system  is  such  that  x 
is  in  the  direction  of  the  flow  and  y  is  in  the  direction  perpendicular  to  x.  The  governing 
equations  for  an  inviscid,  compressible  flow,  in  conservative  form,  can  be  written  as 


au  aE  ar 

dt  ^  dx  ^  dy 


(1) 


(  p  \ 

(  pu  \ 

f  \ 

where  U  =  I 

1  pv 

,E  = 

puv 

\puet  +  ^pu  y 

,F  = 

puv 

p^^  +  ^ 

\pei  / 

\pvet  -f-  “Pu/ 

and  p  =  pT,  et  =  ^T+  +  i;2). 

In  the  above  equations,  u  and  v  are  the  velocities  in  the  streamwise  (i)  and  cross¬ 
stream  (y)  directions  respectively,  p  is  the  density,  T  is  the  temperature,  p  is  the  pressure, 
7  is  the  specific  heats  ratio  and  M\  is  the  upper  stream  Mach  number.  The  equations  have 
been  nondimensionalized  with  respect  to  the  mean  velocity,  temperature  and  pressure  of 
the  upper  stream.  The  height  of  the  channel  is  1.  The  flat  walls  are  located  at  y  =  —0.5  and 
y  =  0.5,  respectively.  Boundary  conditions  at  the  walls  are  that  v  and  normal  derivatives 
of  u,  r  and  p  are  zero. 

For  the  purpose  of  this  paper,  an  explicit  second-order  MacCormack  scheme’"*  is  used 
to  solve  the  above  nonlinear  Euler  equations.  Numerical  solutions  are  advanced  from  time 
f  to  f  -f-  At  by  the  two-step  predictor-corrector  type  algorithm. 


Numerical  simulation  starts  with  an  initially  parallel  flow,  namely  : 

u(x,y)  =  U{y)  =  0.5(1  U2  +  {I  -  U’2)tanh(^)) 
i;{x,y)  =  0 

T(x,y)  =  f(y)  =  -  (/)([/  -  [/jA 

p(x,y)  =  p(y)  = 

T{y) 


(2) 

(3) 

(4) 

(5) 


where  is  the  vOi-ticity  thickness  of  the  mean  velocity  lJ{y)  .  The  temperature  T(y)  is 
obtained  by  using  the  Corroco’s  relation  across  the  stream.  U2  and  T2  are  the  velocity  and 
temperature  of  the  lower  stream,  respectively. 
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To  initiate  the  disturbances,  perturbations  of  u,v,T  and  p  are  superimposed  on  the 
mean  flow  at  the  inflow.  For  the  computations  reported  here,  perturbations  in  the  form  of 
eigenfunctions  were  forced.  Thus  the  boundary  values  at  x  =  0  are  : 


/u(0,y,t)\ 

/Uiy)\ 

/u„(y,w)\ 

v{0,y,t) 

0 

N 

Vn{y,<^) 

7’(0,y,t) 

= 

fiy) 

+ 

fniy,o^) 

p(0,y,t) 
Vp(0,y,t)  / 

p{y) 

\  1  / 

n=I 

\Pn(y,u;)/ 

where  Un  ■,  ,  Tn  ,  Pn  and  p„  are  the  normalized  eigenfunctions  of  the  wave  modes.  The 

normalizations  are  such  that  lun(y,a))|mai  =  1-  N  is  the  number  of  modes  used  in  the 
forcing,  e„  is  the  forcing  amplitude  and  a;  is  the  forcing  frequency.  Computations  of  the 
eigenfunctions  will  be  briefly  discussed  below. 


B.  Eigenfuction  expansion 


For  small  perturbations  forced  at  the  inflow,  the  initial  spatial  growth  is  linear  and 
is  governed  by  the  linearized  Euler  equations.  By  invoking  local  parallelism  of  the  mean 
flow  and  substituting  perturbations  of  the  form  for  u  and  similar  forms  for 

other  variables  into  the  linearized  version  of  (1),  an  eigenvalue  problem  for  the  normal 
modes  is  derived.  In  particular,  it  is  found  convenient  to  form  the  eigenvalue  problem  in 
the  pressure  perturbation.  It  is  straightforward  to  determine  that  the  equation  for  the 
pressure  eigenfunction  is 


cPp  {  2k  dU  ldf\dp  fju-kUyMf  .zV-n 
dy"^  ^  —  ktJ  dy  ^  T  dy  )  dy  ^  \  T  )  ^ 


(7) 


with  boundary  conditions  ^  =  0  at  y  =  —0.5  and  y  =  0.5  For  the  spatial  problem 
considered  here  the  frequency  u>  is  given  and  the  wave  number  k  is  sought  from  (7)  as 
the  eigenvalue.  For  a  closed  domain,  the  eigenvalues  form  a  discrete  set  for  any  given  u;. 
A  typical  dispersion  relation  of  the  eigenvalues  is  shown  in  Figure  2.  The  corresponding 
u{y),  v{y),  T{y)  and  p{y)  can  be  easily  obtained  from  p{y)  by  using  the  linearized  Euler 
equations. 


It  is  expected  that  the  normal  modes  form  a  complete  spectrum  of  small  amplitudes 
waves.  Indeed,  for  bounded  flows,  it  was  shown  that  eigenfunctions  of  hydrodynamic 
instability  problems  form  a  complete  basis  in  the  functional  space'®.  This  fact  have  been 
exploited  in  several  recent  studies  where  eigenfuction  expansions  are  used  in  improving 
numerical  calculations'®  or  in  representing  nonlinear  turbulent  flows'^.  Our  aim  here  is 
to  reconstruct  the  flow  field  obtained  from  the  results  of  direct  numerical  simulations  as  a 
summation  of  eigenmodes.  By  doing  so,  the  spatial  development  of  each  individual  mode 
imposed  at  the  inflow  will  be  established. 
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To  implement  the  expansion  in  eigenmodes,  we  first  record  the  time  series  on  a  mesh 
of  selected  points  For  convenience  of  discussion,  let  the  time  series  of  the  normal 

velocity  at  (xi,yj)  be  j/j, t„).  Here  the  prime  denotes  the  fluctuation  over  its  mean 

value  and  is  the  discretized  time.  Using  a  discrete  FFT,  the  time  series  is  transformed 
into  the  frequency  spectrum,  v*(x,,  t/j,u;).  We  expect  £>'{ j,, , uj)  to  be  a  summation  of 
eigenfunctions  for  fixed  x,  and  u>,  i.e. 

N 

v(xi,yj,u!)  =  '^a„v„(yj,oj)  (8) 

n=  ' 

where  Vn(y,ijj)  are  the  eigenfunctions  at  frequency  w.  The  coefficients  a„  in  (8)  are  found 
by  requiring  that  the  quantity 

J  n 

be  minimum.  That  is,  the  calculated  time  series  spectrum  as  a  function  of  y  is  fitted  in 
the  least-square  sense  by  a  summation  of  normal  mode  eigenfunctions.  Clearly  the  {a„} 
represents  the  wave  amplitude  of  the  corresponding  mode  and  varies  with  x.  It  gives 
a  measure  of  the  strength  of  the  wave  as  it  propagates.  For  the  results  reported  here, 
eigenfunction  expansions  are  performed  at  the  forcing  frequency. 


III.  RESULTS  AND  DISCUSSIONS 

To  demonstrate  the  effects  of  multiple  propagating  wave  modes  and  the  method  of 
eigenfunction  expansion,  we  study  two  cases  here.  In  the  first  case,  a  Class  B  unstable 
wave  and  a  Class  C  neutrally  stable  wave  are  forced  at  the  inflow  boundary.  In  the  second 
case,  two  neutrally  stable  waves,  the  C  and  D  modes,  are  forced.  In  all  the  calculations 
shown  here  a  101  x  201  uniform  mesh  have  been  used  for  the  MacCormack  scheme.  The 
computational  domain  extends  over  8  channel  heights  in  x.  The  vorticity  thickness  of  the 
initial  flow  is  0.02.  The  Mach  number  of  the  ui>per  and  lower  stream  are  Afj  =  3.5  and 
M2  =  1.2  respectively.  The  lower  stream  has  initial  velocity  U2  =  0.29  and  temperature 
T2  =  0.91.  The  two  streams  have  a  speeds  of  sound  ratio  01/02  =  1.2.  The  specific  heats 
ratio  7  =  1.4  for  both  streams.  The  dispersion  relations  of  the  normal  modes  are  those 
given  in  Figure  2. 

A.  CASE  I;  simulation  of  Bi  and  C2  modes 

The  forcing  modes  chosen  here  are  the  Bi  and  C2  modes.  The  forcing  amplitude  is 
0.015  and  forcing  frequency  is  w  =  1  for  each  mode.  From  linear  stability  calculations,  the 
B\  mode  has  a  wave  number  of  2.13  and  spatial  growth  rate  of  0.313.  The  C2  mode  is  a 


neutrally  stable  wave  with  wave  number  3.93.  The  eigenfunctions  of  the  two  inodes  are 
shown  in  Figure  3  where  the  velocity  components  u  and  v  are  pIotte<l  as  functions  c;f  y. 

Figure  4  shows  the  time  series  and  frequency  spectrum  of  the  streamwise  velocity 
fluctuations  at  selected  centerline  locations.  At  the  chosen  forcing  level,  nonlinear  effects 
are  not  developed  within  the  computational  domain.  An  FFT  was  performed  for  the  tim«' 
series  after  t  =  40  when  the  initial  perturbation  has  been  “washed”  out  of  the  domain. 
The  spectrum  shows  a  clear  peak  at  the  forcing  frequency,  u;  =  1.  In  Figure  5  plots  of 
the  root-mean-square  variations  of  velocity  fluctuations  across  the  stream  are  given  for 
selected  downstream  locations.  It  is  seen  that  the  streamwise  velocity  has  a  peak  in  the 
centerline  where  the  shear  layer  is  formed.  However  no  such  peak  is  found  in  the  normal 
velocity  components.  A  comparison  of  r.m.s  fluctuations  with  the  eigenfunctions  given  in 
Figure  3  seems  to  suggest  that  the  n-components  show  more  apparent  modal  differences 
than  the  u-components.  For  this  reason  the  normal  velocity  fluctuations  are  used  for  the 
eigenfunction  expansion  below. 


As  a  measure  of  disturbance  energy,  we  evaluate  the  following  integral  along  the 
downstream  direction, 


e 


(9) 


where  Urms  Kms  r.m.s.  variations  over  the  mean  values.  In  Figure  6{a),  this 

energy  integral  is  plotted  as  a  function  of  downstream  distance  x.  It  is  seen  that  initially 
the  energy  decreases.  Further  downstream  oscillatory  growth  is  found  with  a  period  of 
about  3,3  channel  heights.  We  point  out  that  plots  of  maximum  Upm.i  '■’rma  exhibited 
similar  trends.  In  Figure  6(b),  the  amplitudes  of  B\  and  C2  wave  modes,  obtained  using 
the  least-square  eigenfunction  expansion  described  in  section  II,  are  shown.  (For  the  results 
given,  time  series  at  51  evenly  spaced  y  locations  across  the  stream  have  been  used  for  the 
least-square  fitting.)  It  is  found  that  the  growth  of  B\  wave  over  the  downstream  distance 
fits  extremely  well  with  a  line  of  linear  growth  rate  0.313.  The  amplitude  of  the  C2  wave  is 
shown  staying  more  or  less  the  same  in  the  downstream  direction.  The  slight  variations  of 
the  C2  amplitude  are  due  to  the  accuracy  limitation  of  the  least  square  fitting.  Figure  6(b) 
indicates  that,  although  the  disturbance  energy  in  r.m.s.  measurements  shows  oscillatory 
growth,  each  wave  initiated  upstream  is  still  propagating  linearly.  Clearly  the  co-existence 
of  the  two  waves  leads  to  periodic  cancellation  and  reinforcement  which  results  in  the 
appearance  of  oscillatory  growth  of  the  energy.  Similar  calculations  for  class  A  waves  and 
neutral  waves  have  also  been  carried  out.  However,  since  the  A  waves  have  larger  growth 
rates,  the  oscillatory  behavior  is  not  as  pronounced  as  that  shown  in  Figure  6. 


As  pointed  out  in  the  Introduction,  multiple  wave  modes  are  excited  unless  the  distur¬ 
bance  at  the  inflow  is  exactly  that  of  the  eigenfunction  of  the  chosen  wave  mode.  The  later 
condition  is  usually  difficult  to  meet  in  experiments  and  sometimes  not  always  imposed  in 
numerical  simulations.  For  example,  Lu  &  Wu"*  used  only  the  pressure  perturbation  in  the 
inflow  in  their  simulations.  Often  in  experiments,  the  initial  disturbances  are  excited  by 
random  fluctuations'*  or  “y-component-producing”  element'*.  In  those  studies,  however. 
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the  growth  rate  of  the  unstable  wave  is  relatively  large.  Beeau.se  of  this,  the  rompotu‘nt 
of  the  unstable  wave  grows  quickly  anfl  becomes  dominant  over  other  wave  components 
within  a  short  distance.  As  a  result,  a  linear  growth  rate  in  disturbance  entTgy  is  observed. 
However,  as  shown  in  the  present  calculations,  wlu'n  the  growth  r<ite  of  tin’  linear  insta¬ 
bility  wave  is  relatively  small,  tlu;  other  excited  waves  can  play  a  part  in  the  flow  field  and 
lead  to  an  appearance  of  oscillatory  growth  of  the  disturbance  eiu'rgy. 


B.  CASE  II:  simulation  of  C2  and  modes 

Here  we  study  the  case  when  two  neutral  modes  are  excited.  The  forcing  mod<'.s  are 
the  C2  and  D3  modes.  The  forcing  frecpiency  lj  —  0.25.  The  wave  numbers  of  the  Cb  and 
D3  modes  are  2.76  and  —2.53,  respectively.  The  forcing  amj)litiul('s  are  0.02  for  the  C2 
mode,  and  0.01  for  the  D3  mode. 

In  Figure  7(a)  the  resulting  energy  integral  (9)  as  a  function  of  downstream  distance 
X  is  given.  Using  the  least-square  eigenfunction  exi>ansion,  the  amplitude  of  each  wave  is 
calculated  and  given  in  Figure  7(b).  Figure  7  demonstrates  that,  although  the  energy  of 
the  perturbation  field  has  a  spatial  periodicity,  the  amplitud('  of  each  wave  still  remains 
the  same  as  they  propagate  downstream.  Again  this  is  not  too  surprising  because  the  two 
waves  are  traveling  at  different  phase  velocities  and  amplitude  modulation  is  expected. 
However  it  is  important  to  note  that  the  periodic  structure  of  disturbance  energy  variation 
is  time  independent.  This  point  is  made  clearer  in  figure  8  where  th('  r.m.s.  variation  across 
the  stream  for  u'  and  v'  are  plotted  along  the  downstream  distance.  It  show’s  that,  w’hen 
averaged  in  time,  the  co-existence  of  the  C2  and  D3  waves  product's  a  spatially  periodic 
and  stationary  structure. 

The  period  of  the  structure  can  be  estimated  as  below.  Without  the  loss  of  generality, 
let  two  waves  be  C  cos(kcx  —  wf )  and  D  cos(fcDX  —  cut),  where  kc,  kp  represent  the  w'ave 
numbers  and  C,  D  are  the  amplitudes.  Then  the  mean  square  variation  produced  by  the 
two  waves  is  (overline  denotes  time  averaging)  : 


[C  cos(kcx  —  cut)  -f  Dcos{kQX  —  cut)]2 

=  cos^{kc3^  —  w>t)  -f  2C  D  cos{kc^  —  cut)cos(A:ox  —  cut)  -|-  cos^(fcoj'  —  u.<t) 

=  +  CDcos[(kc  -  kD)x] 

It  is  seen  that  a  streamwise  modulation  is  resulted  with  a  spatial  period  of  Yc-ko  ' 
expression  gives  a  period  of  1.19  for  the  present  case,  which  is  the  same  as  observed  in 
Figure  7  and  8.  In  Ca»se  I,  the  expression  gives  a  period  of  3.49  by  substituting  D  by  B\ 
in  (10). 
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IV.  CONCLUSION 


One  differenee  betwe>cn  the  numerical  simulation  of  bounded  and  uul)ounded  slusir 
layers  is  tliat  the  former  permits  multiple  propagating  wave  modes.  This  becomes  iinpoi' 
tant  when  the  growth  rates  of  the  unstable  Wcives  are  small.  It  was  showm  here  that  two 
forced  modes  can  produce  oscillatory  growth  in  disturbance  energy  even  though  each  indi- 
vid\ial  mode  is  propagating  linearly.  In  particular,  it  was  shown  that  two  excited  neutral 
modes  with  different  phase  velocities  can  lead  to  a  station:iry  and  periodic  structure  cT  th(' 
r.m.s  fi\ictuations.  Moreover,  a  method  of  effectivily  detecting  and  separating  th<'  growth 
of  each  individTial  wav'e  mode  is  given.  Finally  it  was  shown  numerically  that,  by  fitting 
the  disturbance  field  as  a  sunmiation  of  eigenfTinctions.  the  amplitude  of  t  ach  mode  can 
l)e  accurately  monitored. 
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Schematic  drawing  of  a  bounded  shear  layer 
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Figure  2.  Dispersion  relation  of  the  normal  modes  associated  with  a  bounded  supersonic 

shear  layer.  - unstable  modes, - neutral  modes.  Mode  number  is  denoted  by  subscripts. 

w  is  the  frequency,  and  k^  are  the  real  and  imaginary  parts  of  the  wave  number,  respectively. 
M,  =  3.5,  Mi  -  1.2,  a,/oj  =  1.2. 
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x  =  6.4 


^’rms  x-O.O  x=1.2  x=3.2 


Figure  5.  Root-mean-square  variations  across  the  stream  at  indicated  downstream  locations 
The  velocity  components  are  shown.  B\  and  Cj  modes  are  forced  at  the  inflow. 


X 


Figure  6.  (a)  Energy  integral  (9)  as  a  function  of  downstream  distance  x,  showing  oscillatory 

growth.  Vertical  dotted  lines  have  separation  distance  of  3.3.  (b)  Amplitudes  of  Bi  and  Ci  modes 
calculated  using  the  least-square  fitting  with  eigen funcions.  The  solid  straight  line  represents  an 
exponential  growth  rate  of  0.313.  ,  p 


1  2  3  4  5  6  7  8 
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Figure  7.  (a)  Energy  integral  (9)  as  a  function  of  downstream  distance  z.  The  C2  and  -D3  modes 

were  forced.  Vertical  dotted  lines  have  separation  distance  of  1.19.  (b)  Amplitudes  of  the  Cj  and 

Dz  modes  calculated  using  least-square  fitting  with  eigenfunctions. 
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Figure  8.  Root-mean-square  fluctuation  as  function  of  »  and  y.  (a)  Streamwise  velocity 
component;  (b)  normal  velocity  component;  (c)  contours  of  the  normal  velocity  component  shown 
in  (b). 
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